arXiv:1502.00904vl [astro-ph.SR] 3 Feb 2015 


Research in Astronomy and Astrophysics manuscript no. 
(DTpX: ms2059.tex; printed on February 4, 2015; 1:23) 


Light curve solutions of six eclipsing binaries at the lower limit of 
periods of the W UMa stars 

Diana P. Kjurkchieva 1 , Dinko P. Dimitrov 2 and Sunay I. Ibryamov 1,2 

1 Department of Physics, Shumen University, 115, Universitetska Str., 1712 Shumen, Bulgaria; 
d. kyurkchieva@ shu-bg. net 

2 Institute of Astronomy and National Astronomical Observatory, Bulgarian Academy of Sciences, 72, 
Tsarigradsko Shose Blvd., 1784 Sofia, Bulgaria 

Received ; accepted 


Abstract Photometric observations in V and I bands of six eclipsing binaries at the lower 
limit of the orbital periods of W UMa stars are presented. Three of them are newly discovered 
eclipsing systems. The light curve solutions revealed that all short-period targets were contact 
or overcontact binaries and added new six binaries to the family of short-period systems with 
estimated parameters. Four binaries have equal in size components and mass ratio near 1. The 
phase variability of the V-I colors of all targets may be explained by lower temperatures of 
their back surfaces than those of their side surfaces. Five systems revealed O’Connell effect 
that was reproduced by cool spots on the side surfaces of their primary components. The 
light curves of V1067 Her in 2011 and 2012 were fitted by diametrically opposite spots. The 
applying of the criteria for subdivision of the W UMa stars to our targets led to ambiguous 
results. 

Key words: methods: data analysis, catalogs, stars: fundamental parameters, stars: bina¬ 
ries: eclipsing: individual: 1SWASP J173828.46+111150.2,1SWASP J174310.98+432709.6 
= V1067 Her, NSVS 11534299, NSVS 10971359, NSVS 11234970, NSVS 11504202 

1 INTRODUCTION 

Most of the W UMa stars consisting of solar-type components have orbital periods within 0.25 d < P < 
0.7 d. They are recognized by continuous brightness variations and nearly equal minima. The short orbital 
periods of these binaries mean small orbits and synchronized rotation and orbital revolution. 
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Table 1 List of our targets 


Target Name 0:2000 82000 Period [d] V [mag] AV Ref 


1 

1SWASP J173828.46+111150.2 17 38 28.459+11 11 50.05 0.249383 

13.80 

0.24 

1 

2 

1SWASPJ174310.98+432709.6 17 43 10.977 +43 27 09.48 0.258108 

13.20 

0.44 

2 

3 

NSVS 11234970 

19 32 07.875 +14 58 26.90 

0.25074 

13.66 

0.63 

3 

4 

NSVS 11504202 

20 27 27.304+12 30 50.66 

0.24614 

13.67 

0.66 

3 

5 

NSVS 11534299 

20 48 20.859 +06 59 22.98 

0.22465 

13.40 

0.21 

4 

6 

NSVS 10971359 

18 29 20.479 +07 58 07.26 

0.27977 

13.31 

0.38 

3 


References: 1 - Lohr et al. (2012); 2 - Blattler & Diethelm (2000); 3 - this paper; 4 - Hoffman et al. (2009) 

The statistics of W UMa stars around the lower limit of periods is quite poor (Terrell et al. 2012) because 
the period distribution of binaries reveals a very sharp decline below 0.27 days (Drake et al. 2014). Another 
reason is that the short-period binaries consist of late stars and thus are faint objects for detailed study. 

However, the short-period contact systems are important objects for the modern astrophysics at least 
in two ways: (a) the empirical relation period-luminosity of W UMa’s allows to use them as distance and 
population tracers; (b) these binaries are natural laboratories for study of late stage of the stellar evolution 
connected with the processes of mass and angular momentum loss, merging or fusion of the stars, etc. 

Lately, the interest to binaries around the low-period limit increases because they provide constraints 
on the theories for the formation and evolution of the low-mass stars. Moreover, one of the hypotheses for 
the origin of the hot Jupiters is that they are products of proto-planetary disk formed by merging of short- 
period, low-mass binaries via ’’magnetic braking” if substantial amount of mass remains in orbit around the 
primary (Martin et al. 201 1). 

Fortunately, the modern large stellar surveys during the last decade allowed to discover binaries with 
shorter and shorter periods Rucinski (2007), Pribulla et al. (2009), Weldrake et al. (2004), Maceroni & 
Rucinski (1997), Dimitrov & Kjurkchieva (2010), Norton et al. (2011), Nefs et al. (2012), Davenport et al. 
(2013), Lohr et al. (2014), Qian et al. (2014), Drake et al. (2014). The majority of the newly discovered bi¬ 
naries from space missions ( Kepler, C 0 R 0 T, etc.) and ground-based projects (ASAS, SuperWASR Catalina, 
LINEAR, NS VS, etc.) were classified as contact systems. But most of them need follow-up observations 
and study. 

In this paper we present photometric observations and light curve solutions of six binaries around the 
lower limit of the orbital periods of W UMa stars. We chose to study such systems because the statistics of W 
UMa stars with periods around quarter day is quite poor (Terrell et al. 2012). Three of our targets are known 
binaries: two binaries (1SWASP J173828.46+111150.2,1SWASP 1174310.98+432709.6= V1067 Her) are 
from the SuperWASP photometric survey (Pollacco et al. 2006) and one system (NSVS 11534299) is from 
the NSVS data base (Wozniak et al. 2004). Three of our targets (NSVS 11234970, NSVS 11504202, NSVS 
10971359) are newly discovered binaries. Table I presents the coordinates of our targets and information 
for their light variability. 
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Table 2 Journal of the Rozhen photometric observations 


Target Date Filter Exposure [sec] Phase range 


1 

2011 July 10 

v,i 

120,120 

0.21-1.17 


2011 Aug 03 

VI 

120,120 

0.67-0.89 


2011 Aug 21 

VI 

120,120 

0.60-0.94 


2012 June 02 

VI 

120,120 

0.61-1.41 


2012 June 03 

VI 

120,120 

0.24-1.35 


2012 July 07 

VI 

120,120 

0.90-1.74 

2 

2011 Aug 11 

VI 

120,120 

0.20-1.08 


2011 Aug 19 

VI 

120,120 

0.83-1.74 


2012 June 04 

VI 

120,120 

0.86-1.33 


2012 July 08 

V,I 

120,120 

0.35-1.38 

3 

2012 Aug 04 

VI 

90,120 

0.60-1.94 

4 

2012 Aug 05 

VI 

120,120 

0.63-1.91 

5 

2012 Aug 06 

VI 

120,120 

0.19-1.32 

6 

2011 July 08 

V,I 

120,120 

0.45-1.24 


2011 July 11 

V,I 

120,120 

0.16-0.90 


2 OBSERVATIONS 

The CCD photometry of the targets was carried out in 2011—2012 at Rozhen National Astronomical 
Observatory (Bulgaria) with the 60-cm Cassegrain telescope using the FLI PL09000 CCD camera (3056 x 
3056 pixels, 12 fj,m/ pixel, field of 17.1'x 17.1'). The average photometric precision per data point is 0.005 
mag in I band and 0.008 mag in V band. Table 2 presents the journal of our simultaneous VI observations. 
In addition we obtained several observations in BVI colors of the targets 1, 3 and 4. 

The photometric data were reduced by IDL package (subroutine DAOPHOT). We used more than three 
standard stars in the observed fields (Fig. 1). Table 3 presents their colors: I from USNO-Bl.O catalogue 
and V from GSC 2.3.2 catalogue. The magnitudes of the targets in Table 3 correspond to their out-of-eclipse 
levels from our observations. 

We carried out preliminary time-series analysis of our data by the software PerSea. The comparison of 
the results with the published ones (Table 1) revealed two discrepancies. 

(a) The periods of targets 1 and 2 phased well our data. But the period 0.24914 d of target 5 (see 
further Table 5) we obtained was around 11 % longer than the previous value (Table 1). This result is 
not so surprising because the follow-up observations of variable stars from different data bases sometimes 
reveal different periods than the previous values (Norton et al. 2011; Mighell & Plavchan 2013; Lohr et al. 
2013; Kjurkchieva et al. 2014; Mayangsari et al. 2014). This is a consequence of the automated frequency 
analyses of the huge data sets. They based on different approaches which turned out with different precision 
for different types of variabilities. Moreover, sometimes they give several periods for each target (as well as 
false alarm detections of exoplanets). 
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Fig-1 The fields of the targets 

(b) The amplitudes of variability of targets 1 and 2 turned out considerably bigger than the previous 
ones. We assume that the reason for this discrepancy is the low space resolution of the previous photomet¬ 
ric observations by small telescopes (13.7"pixel _1 for SWASP survey and 14.4"pixel _1 for NSVS) that 
prevents separation of two neighboring stars. 
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Table 3 Magnitudes of the targets and their standard stars 

Star GSC-ID I V B 

(USNO B1) [mag] [mag] [mag] 

targ 1 


stl 0100100125 
st2 0100101007 
st3 0099702419 
st4 0099702343 
st5 0099702387 
targ 2 

stl 0310001679 
st2 0310001797 
st3 0310001604 
targ 3 

stl 1049-0454686 
st2 1049-0454529 
st3 1049-0454842 
targ 4 

stl 0109500807 
st2 0109501223 
st3 0109500601 
st4 0109501589 
st5 0109500781 
targ 5 

stl 0052402305 
st2 0052401323 
st3 0052401783 
targ 6 

stl 0102300733 
st2 0102302368 
st3 0979-0436284 


12.21 

13.46 

14.68 

11.55 

12.67 

13.41 

10.83 

12.20 

13.62 

11.04 

12.49 

14.15 

11.10 

12.57 

14.09 

11.85 

13.12 

14.10 

12.04 

12.82 


12.49 

13.00 


13.45 

13.66 


11.13 

12.01 


12.45 

13.66 

14.33 

12.88 

13.44 

13.87 

10.91 

12.45 

13.33 

11.75 

12.88 

13.41 

12.43 

13.67 

14.64 

12.63 

13.19 

14.79 

10.98 

12.27 

13.70 

11.22 

12.78 

14.65 

10.66 

11.33 

11.93 

10.88 

12.38 

14.02 

11.97 

13.19 


10.71 

12.13 


11.50 

11.93 


11.94 

13.32 


12.92 

13.31 


12.14 

12.34 


10.87 

11.93 


13.77 

13.97 



The close nonvariable star may enter the same pixel or aperture during the photometric measurements 
of the variable and thus to reduce its amplitude of variability. 

3 LIGHT CURVE SOLUTIONS 

The shapes of the light curves (Figs. 2, 3, 4, 5, 6, 7) implied that our targets are nearly contact or overcontact 
systems that was expected for their short orbital periods. The big amplitudes of their light variabilities mean 
that they are caused by eclipses. 

It is well known that the determination of the mass ratio through light-curve solution is an ambiguous 
approach compared with that by radial velocity solution. 
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Table 4 The 2 MASS color indices J-K and corresponding mean 
temperatures T m of the targets 


Target 1 2 3 4 5 6 

J-K 0.645±0.033 0.590±0.024 0.749±0.030 0.750±0.026 0.558±0.034 0.618±0.031 
T m 4720±120 4930±140 4320±130 4320±110 5090±220 4810±130 


However, the rapid rotation of the components of the short-period binaries is serious obstacle to obtain 
precise spectral mass ratio from measurement of their highly broadened and blended spectral lines (Bilir 
et al. 2005, Dali & Schmidtobreick 2005). On the other hand, their eclipse depths depend strongly on the 
potentials and the mass ratios. 

Taking into account these considerations we solved the Rozhen light curves of the six targets using the 
code PHOEBE (Prsa & Zwitter 2005) by the following procedure. 

We determined in advance the mean temperatures T m of the binaries (Table 4) by their infrared color 
indices (J-K) from the 2MASS catalog and the calibration color-temperature of Tokunaga (2000). 

At the first stage we adopted Ti=T m (a good approximation for contact binaries with close temperatures 
of the components) and searched for solutions for fixed 7) varying the initial epoch To, period P, secondary 
temperature X 2 , orbital inclination i, mass ratio q and potentials O] . 2 - 

We adopted coefficients of gravity brightening 31 = <72 = 0.32 and reflection effect A\ = A 2 = 0.5 
appropriate for late stars while the limb-darkening coefficients for each component and each color were 
auto-updated to the stellar temperature according to the tables of Van Hamme (1993). 

In order to reproduce the O’Connell effect we added cool spots on the stellar surfaces and varied spot 
parameters (longitude A, latitude /3, angular size a and temperature factor k). Moreover, the best fit in the 
two colors required a small contribution of third light ( 3 . 

Finally, we searched for the best fits for fixed q, spot parameters and third light contributions by adjusting 
the primary temperature 7 j slightly above T m and changed correspondingly 7 2 , V. [ .2 and i. 

Table 5 contains the results of our light curve solutions. The parameter errors are the formal 
PHOEBE errors. The synthetic light curves corresponding to the parameters from Table 5 are shown in 
Figs. 2, 3, 4, 5, 6 , 7, 8 as continuous lines. 

The results of our light curve solutions can be used for estimation of the global parameters of the 
targets (masses, radii and luminosities) using statistical relations and supposition that their components are 
approximate MS stars. 

4 ANALYSIS OF THE RESULTS 


The analysis of the light curve solutions led to several conclusions. 
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Fig. 2 Light curves of target 1 in V and I bands (the V data are shifted vertically by l m for better 
visibility) and their fits 



Fig. 3 Light curves of target 2 in V and I bands in 2011 (the V data are shifted vertically by 
0.3 m for better visibility) and their fits 

(1) All targets are almost contact (~CB) or overcontact binaries (OC). This result is expected taking 
into account the short periods of the binaries (Table 1). 

(2) The q values are near 0.5 for two binaries (targets 1 and 6) and near 1 for the rest ones. This result 
refers to our sample only and should be assumed as accidental because the q values of W UMa stars are 
rather evenly distributed (fig. 1 in Csizmadia & Klagyivik 2004). 

(3) The stellar components are of K spectral type. This is a natural consequence of the average infrared 
colors of the targets (Table 4). 

(4) The temperature differences of the components of the binaries are below 630 K. The small tempera¬ 
ture differences are expected for contact or overcontact systems and means that their components are almost 
in thermal contact. 

(5) Four targets (2, 3,4, 5) have equal in size components. Just they have mass ratio near 1. Hence, their 
components obey the same mass-radius relation. 
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Fig. 4 Light curves of target 2 in V and I bands in 2012 (the V data are shifted vertically by 
0.3 m for better visibility) and their fits 



Fig. 5 Light curves of target 3 in V and I bands (the V data are shifted vertically by l m for better 
visibility) and their fits 

(6) Five targets revealed O’Connell effect that was reproduced by cool spots (Table 6) on the side sur¬ 
faces of their primary components. The light curves of target 2 in 2011 and 2012 were fitted by diametrically 
opposite spots (visible at the second and first quadratures correspondingly). This result can be explained by 
differential rotation or by new spot cycle. But the equal sizes and temperatures of the two spots (Table 6) 
support rather the first supposition. 

The modeling of the O’Connell effect of late stars by cool spots is preferred approach of almost all 
authors. Cloud of circumstellar, absorbing matter that orbits the corresponding star may reproduce the 
O’Connell effect as the cool photospheric spot. Only photometric data (as the ordinary case is) cannot 
resolve this ambiguity and one adopts modeling by spots. 

(7) The cool spot of target 5 was invoked to reproduce the distortions of its light curve at the minima 


and maxima. 
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Fig. 6 Light curves of target 4 in V and I bands (the V data are shifted vertically by l m for better 
visibility) and their fits 



Fig. 7 Light curves of target 5 in V and I bands (the V data are shifted vertically by l m for better 
visibility) and their fits 

(8) The V-I indices of the targets undergo orbital variability with reddening at the light minima (Fig. 9). 
The amplitudes of the V-I curves are proportional to the light amplitudes (Table 7). This result means that 
the temperature of the back surface of the close stellar components is lower than those of their side surfaces. 

(9) No target of our sample was identified as an X-ray source and H (1 emission source. Thus the low 
activity of the W UMa stars was confirmed. 

5 SUBCLASSIFICATION OF OUR TARGETS 

The contact binary stars were divided into two sub-types, A and W, according to the following criteria: (a) 
the ratio RJT (radius to temperature): the larger star is the hotter one for A-type system while the smaller 
star is the hotter one for W-type system (Binnendijk 1970); (b) temperature T or spectral type: the A type 
systems are earlier than the W type binaries which components are of G and K spectral type; (c) period P: 
the W-type binaries have shorter periods of 0.22 to 0.4 day (Smith 1984); (d) mass ratio q\ values of A type 
are smaller than those of W type. 
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Fig. 8 Light curves of target 6 in V and / bands (the V data are shifted vertically by 0.2 m for 
better visibility) and their fits 


Table 5 Parameters of the light curve solutions of the targets 


Parameter 

1 

2 

3 

4 

5 

6 

T 0 2450000+ 

6091.94682 

5785.294939 

6144.364434 

6145.36124 

6146.46937 

5751.49320 


±0.00021 

±0.000151 

±0.000208 

±0.00011 

±0.00021 

±0.00021 

P 

0.249383 

0.258108 

0.25074 

0.24614 

0.24914 

0.27977 




±0.00027 

±0.00024 

±0.00057 

±0.00041 

q 

(N 

O 

O 

d 

-H 

00 

00 

d 

1.009±0.005 

0.986±0.006 

VO 

O 

O 

d 

-H 

0 

00 

05 

d 

0.872±0.003 

0.462±0.002 

i 

66.85±0.09 

75.68±0.21 

74.36±0.26 

79.4±0.2 

56.3±0.4 

66.8±0.2 

Ti 

4900±48 

5152±94 

4600±90 

4402±55 

5150±61 

4936±48 

t 2 

4544±28 

4626±56 

3969±70 

4219±54 

5014±45 

4423±37 

Ql 

2.843±0.004 

3.640±0.012 

3.615±0.009 

3.717±0.013 

3.548±0.014 

IT) 

O 

O 

d 

-H 

to 

00 

(N 

^2 

2.843±0.004 

3.640±0.012 

3.615±0.009 

3.717±0.013 

3.550±0.008 

10 

0 

0 

d 

-H 

to 

00 

(N 

ri 

0.445 

0.400 

0.399 

0.363 

0.389 

0.451 

T2 

0.319 

0.399 

0.396 

0.338 

0.365 

0.317 

h 

0.746 

0.643 

0.715 

0.598 

0.566 

0.787 

h 

0.254 

0.357 

0.285 

0.402 

0.434 

0.213 

I 2 /I 1 

0.340 

0.555 

0.399 

0.672 

0.767 

0.271 

fillout 

0.033 

0.229 

0.209 

-0.045 

-0.0024 

0.061 

fillout 

0.033 

0.229 

0.209 

-0.045 

-0.003 

0.061 

configuration OC 

OC 

OC 

~CB 

~CB 

OC 


Table 6 Parameters of the cool spots on the targets 


Target 

0 

1 

122±1 

22011 

90±1 

2 2 012 

90±1 

3 

90±1 

4 

90±1 

5 

40±1 

6 

90±1 


\ a k I3 

297±1 17±1 0.83±0.01 0.029(1) 
90±1 20±1 0.80±0.01 0.066(1) 
270±1 20±1 0.80±0.01 0.034 (I) 
270±1 18±1 0.90±0.01 0.055 (I) 
270±1 10±1 0.90±0.01 0.085 (I) 
Oil 20±1 0.88±0.01 0.073 (V) 
90±1 15±1 0.90i0.01 0.057(1) 
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Phase 


Fig. 9 V — I light curve of target 1 

Table 7 Amplitudes of the reddened (V — I) colors and ampli¬ 
tudes of I variability of the targets 

Target 1 2 3 4 5 6 

A(V - I) 0.06 0.08 0.10 0.14 0.03 0.07 
A / 0.30 0.55 0.55 0.55 0.21 0.32 


Table 8 Subclassification of the targets 


Criterion/target 1 

2 

3 

4 

5 

6 

T/R 

A 

A 

A 

A 

A 

A 

T 

W 

w 

W 

W 

W 

W 

P 

W 

w 

W 

W 

W 

W 

9 

A 

w 

w 

w 

w 

A 

9, h/h 

a : 

H/A 

7 

H/A 

H/W A 


Csizmadia & Klagyivik (2004) introduced H subtype systems (H/A and H/W) with large mass ratio 
(q > 0.72) which energy transfer is less efficient than that in other types of contact binary stars. They found 
that the different subtypes of W UMa’s are located into different regions on the mass ratio - luminosity 
ratio diagram but above the line A = q 4A> (A = fo/h) representing the mass-luminosity relation for the MS 
stars (fig. 1 in Csizmadia & Klagyivik 2004). 

The classification of our targets according to the foregoing criteria is presented in Table 8. It reveals that 
almost all targets belong simultaneously at least to the two subclasses. Only target 3 is deeply below the 
line representing the mass-luminosity relation for the MS stars (uninhabited place) This result means that 
the proposed criteria and subdivision are ambiguous and, of course, that the stellar world is richer than one 
expects. 
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6 CONCLUSION 

We obtained light curve solutions of six binaries with periods around 0.25 d, the lower limit of the periods 
of the W UMa stars. Three of the targets were newly discovered eclipsing systems. 

The solutions revealed that all systems are contact or overcontact binaries. 

The colors of the investigated systems become redder at the light minima. Such a behavior may mean 
that the temperatures of the back surfaces of the close components are lower than those of their side surfaces. 

The mass ratios of the targets are grouped around values 0.5 and 1 and those with mass ratio near 1 have 
equal in size components. 

Five targets revealed O’Connell effect that was reproduced by cool spots on the side surfaces of their 
primary components. The light curves of V1067 Her in 2011 and 2012 were fitted by diametrically opposite 
spots. 

The applying of the criteria for subdivision of the W UMa stars to our targets led to ambiguous results. 
New criteria are necessary for the subdivision of the numerous W UMa systems. 

This study adds new six systems with estimated parameters to the family of short-period binaries. They 
could help to improve the statistical relations between the stellar parameters of the low-mass stars. 
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